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1. Objectives of Research 


This is the final report on 2.5 calendar years of contract research on 
magnetosphere imaging in the extreme ultraviolet. The objectives of the 
research contract can be summarized in three related items: 

Objective 1; Investigate the feasibility of the EUV imaging concept 
by using Lockheed in-situ ion measurement data. 

Objective 2; Investigate the technical resources required to make 
scientific interpretations of such EUV images by end-to-end simulations of 
the imaging concept. 

Objective 3: Construct and test the key element for EUV 
magnetosphere imaging, the multilayer mirror with normal-incidence 
reflectance designed for the selected solar resonance line. 

2. Accomplishments 

Researches as per contract have successfully fulfilled all three 
scientific objectives stated above. The body of this report emphasizes 
accomplishments in the second year. Accomplishments in the first year 
are summarized in the Annual Progress Report for 1991, which is included 
here as Appendix 1. 

Objective 1: Early in the contract period we published a paper on the 
feasibility of magnetosphere imaging in the EUV in the Geophysical 
Research Letters (Chiu et al., 17, 267; see Exhibit 1 in Appendix 1). 
Subsequently, we follow up the work with a more detailed work which has 
been accepted by Planetary and Space Science (see Exhibit B of this 
report). Both papers are now referenced repeatedly by workers in the 
field of EUV imaging of the magnetosphere. Instrument groups are now 
using its results to determine source requirements for instrument design 
proposals for NASA's prospective Small Explorer and Inner Magnetosphere 
Imager missions. 

Objective 2: Unlike the more straight forward estimates of photon- 
counts based on in-situ ion flux data for accomplishing Objective 1, this 
objective is much more difficult. It attempts to answer the question: How 
does one convert morphological features of EUV photon images into ion 
distributions? Here, subtle effects such as Doppler shifts due to thermal 
and non-thermal ion motion in a magnetic field become important and 
must be simulated. Such simulations require considerable computational 
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resources and are accomplished in 1992. The result is a paper recently 
submitted to the Geophysical Research Letters (see Exhibit A of this 
report), since our colleagues felt that the results are timely and important. 
This paper points out that Doppler shift effects in a magnetic field must be 
taken into account in order to convert EUV morphology into ion density 
distribution/morphology. 

Aside from fulfilling Objectives 1 and 2 by publication of three 
science papers, at least five invited papers have been given, by 
investigators of this contract, in American Geophysical Union, International 
Society for Optical Engineering and COSPAR meetings, as well as numerous 
contributed papers. The principal investigator has chaired and organized 
two magnetosphere imaging special sessions for scientific society meetings. 
All such activities were partially under the auspices of this contract. 

Objective 3: To image the magnetosphere with an EUV instrument, it 
is not enough to simply have an EUV normal-incidence multilayer mirror 
tuned to the proper wavelength. The magnetosphere is a large object, so it 
is necessary to demonstrate that a multilayer mirror with a wide field of 
view can be successfully constructed and tested. This is not a simple task 
because it requires proper wavelength band transmission for a highly 
curved surface on a mirror of large diameter. To demonstrate it is feasible, 
we have constructed such a mirror at 304 angstrom (the helium ion line) 
and demonstrated its reflectivity. The mirror constructed under auspices 
of this contract is indeed the world's largest 304 angstrum multilayer 
mirror to date. The detailed report of the construction and testing is 
included below as part of the main body of this report: 
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Fabrication and Test of a Wide Field Multilayer Mirror 


The objective of the laboratory research conducted under this contract was 
to fabricate and test a multilayer deposited on a highly curved mirror surface. 
The multilayer was designed for normal incidence reflection of the emission 
line of He II at a wavelength of 304 A. The mirror was of a design suitable for 
broad field imaging of the reasonantly scattered solar He line from He+ in the 
earth's plasmasphere. The mirror must have a large field of view f approximately 
40 deg full width, in order to effectively image large scale structures in the 
plasmasphere. This requires a short focal length to match available detector 
sizes and to minimize telescope weight and volume. 

A spherical protype mirror was chosen having a radius of curvature of 9.8 
cm (focal length of 4.9 cm) and a diameter of 14 cm. The sagitta for this 
highly curved mirror is 2.8 cm and the angle that the mirror surface makes with 
its axis varies from 90 deg to 45 deg, center to edge. Since the multilayer 
deposition rate is both a function of distance and angle of incidence, this 
poses a challenge to produce a multilayer that is uniform in response over the 
mirror's surface. 

A multilayer consists of interleaved thin films of materials with 
differing optical constants. For normal incidence reflection of 304 A photons, 
molybdenum and silicon were chosen. A vacuum sputtering process, involving 
cylindrical planar magnetrons, is used to deposit thin films of these metals. 
These devices generate a magnetically confined plasma in an argon atmosphere at 
a pressure of -2 mtorr. An applied electric field accelerates ionized Ar atoms 
which strike a target of the material to be deposited. The Ar atoms have 
sufficient energy to knock out, or sputter, atoms of the target material that 
are subsequently deposited on the multilayer substrate. 

The multilayer deposition system is shown in Fig. 1. Two magnetrons, one 
with a Si target and the other with Mo, appear as the the cylinders that are 
diametrically opposed just below the disk near the top of the photograph. The 
deposition regions of magnetrons are bounded by the cylinders and the disk. 
Multilayer substrates may be mounted in holes cut in the disk, whose rotation 
then carries the substrate alternately over the Si and Mo deposition regions to 
fabricate the multilayer. Precise calibration of the deposition rates and 
maintaining them constant in time allows layer thicknesses to be selected and 
controlled throughout the multilayer. Areal uniformity of the multilayer is 
obtained by using appropriately shaped masks between the magnetrons and the 
substrate which maintain constant deposition rates over its entire area. 

In order to determine the required mask shapes for deposition on the 
highly curved mirror, an aluminum model of the mirror was fabricated. This 
allowed measurements to proceed without placing the mirror at risk. Since the 
mirror is symmetrical, it was decided to rotate the substrate about its axis at 
1 rps during the depositions. Six facets, approximately one cm square, were cut 
from a polished silicon wafer and mounted radially along the mirror model's 
surface. The model had been made with slightly larger radius of curvature to 
compensate for the 0.2 mm facet thickness. Thus the surfaces of the facets 
closely simulate the shape of the mirrors surface. As a means of calibrating 
the deposition rates a strip of photo-resist was applied across each facet. 

The faceted model was rotated over a magnetron's deposition region for a 
measured time interval and the photo-resist was subsequently desolved to 
produce a step equal in height to the amount of material deposited. Measurement 
of the step height and the deposition time calibrated the rate. This was done 
for each facet as a function of radius on the mirror model. Iteration of this 
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process with new facets for each trial and magnetron, while varying the shape 
of the mask, was carried out until satisfactorily constant step heights were 
obtained for all facets. The final step heights of approximately 1200 A are 
shown in Fig. 2. The deposition times of 1200s and 600s for the Mo and Si 
magnetrons have been selected to give a measurement accuracy for the step 
heights of approximately +- 1.5 %. The final peak to valley spread for Mo is 
1.4%, while that for Si is 4.8%. Deposition rates calculated from these data 
are used to determine the dwell time over each magnetron that is required to 
fabricate the selected multilayer design. 

Because of spherical aberration in a mirror with such a small focal ratio 

(f/0.35), only a narrow annulus of the entrance aperture can be used if 

acceptable angular resolution is to be obtained over a 40 deg field. While this 
reduces the usable effective area of the mirror, it has a beneficial effect, 
since baffling the mirror in this manner restricts the angle of incidence onto 
the mirror to a very narrow range. One may then tailor the Bragg peak of the 
multilayer to the average angle and the narrow angular range provides more 
uniform reflectivity over the field of view. A ray trace of the optics to 
determine this average angle of incidence was beyond the scope of this 

investigation. However, parameters of the present mirror are very close to 

those of telescopes on the ALEXIS instrument (Priedhorsky et al . , 1988 and 
Bloch et al., 1992). In that case the average incident angle was 23 deg with a 
spread of ~ +— 3 deg and these parameters were used in this investigation. 

The response of the selected multilayer design is shown in Fig. 3. This 
plot is the Bragg reflectivity of the multilayer for 304 A photons as 
calculated from optical constants of Mo and Si (Windt, 1989) . The multilayer 
has a 2d-spacing of 368 A and 5 layer-pairs with Mo being 36% of the layer-pair 
thickness. In predicting the 2d spacing on each facet from the data of Fig. 2, 
a peak to valley spread of 2.9% was calculated. This is less than the 
percentage spread in Si step heights, 3ince the lowest rate in Si is partially 
compensated by the highest rate in Mo at a radius of 1.25 cm and thus reduces 
the maximum spread. At the same magnetron deposition rates determined from data 
of Fig. 2, dwell times of 55 s over the Si magnetron and 67 s over that of Mo 
were determined. 

Multilayers were fabricated on silicon wafer facets located, as before, at 
the six radii on the model mirror used for the rate calibrations. The measured 
304 A reflectivity of each facet as a function of angle of incidence is shown 
in Fig. 4. There are clearly non uniformities from facet to facet. We attribute 
this to the differences in each facet's exposure as it enters and leaves the 
deposition region. These differences were negligable in the 600 s and 1200 s 
deposition times used in the rate calibrations but are apparent in the shorter 
exposures used in the multilayer fabrication. However, the least variation 
occurs near the angle of 23 deg. which is the average incident angle a 
telescope using this mirror would operate. Also, in a working telescope, the 
center of the mirror is occulted by the detector and is thus not illuminated. 
The variation of response is smaller if only the outer three facets are 
considered. Angular responses for these three facets are shown in Fig. 5 and 
data on average reflectivities for a +- 3 deg angular spread are given in Table 
I for both the full aperture and neglecting the three central facets. Maximum 
reflectivity variation for full aperture of ^ M occurs at 26 deg. For the 
limited aperture the maximum is h — 4 % at 20 deg with a variation of less than 
+- 1 % and +- 3 % at the other angles. 


TABLE I: AVERAGE REFLECTIVITY +- HALF OF PEAK TO VALLEY SPREAD 


All Facets 
Outer Facets 


20 

.146+-. 007 
. 143+-. 006 


Angle (deg) 
23 

.146+-. 005 
. 148+-. 001 


26 

.138+-. 010 
.142+-. 004 


The data of Fig. 4 have been fitted by curves similar to Fig. 3 by varying 
the 2d-spacing and Mo fraction. The fits are shown in Fig. 6 and the results 
summarized in Table II. 


TABLE II: DERIVED 2D-SPACINGS AND MO FRACTIONS FOR TEST FACETS 


Radius 

2d 

Mo Fraction 

(cm) 

(A) 


0.0 

358 

.36 

1.25 

365 

.40 

2.5 

367 

.33 

3.75 

371 

.33 

5.0 

376 

.32 

6.25 

369 

.36 

>le II show 

a peak 

to valley s; 


Mo fraction. This large variation in Mo fraction is much greater than expected 
and is not presently understood. It is unlikely that the physical thickness of 
Mo varies by this much and the Si layer thickness varies inversely to 
compensate, so the 2d-spacing (layer-pair thickness) is constant to within the 
observed 5 %. 


The variation in reflectivity indicated in the limited aperture data of 
Table I is sufficiently small that it was decided to deposit a multilayer on 
the mirror using the same procedures as for the model and facets. This 
deposition was completed and a photograph of the multilayer mirror is shown in 
Fig. 7. 

Because of the short focal length and small radius of curvature of the 
mirror, it is very difficult to test. A system was set up using a 
micro-channel plate to measure the reflectivity as a function of wavelength at 
a single location on the mirror. These measurements, along with similar 
measurements for the multilayer facet at a radius of 6.25 cm, are shown in Fig. 
8. The response calculated from optical constants and the parameters of Table 
II for the r=6 . 25 cm facet is also shown by the solid line in the plot. There 
is reasonable agreement at the longer wavelengths but the mirror data point at 
304 A and the facet point at 256 A lie off the curve outside of counting 
statistics, which are approximately the the size of the symbols. 
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In summary, we have successfully demonstrated the ability to deposit a 
multilayer, with uniformity better than 4 %, on highly curved mirror surfaces 
that are suitable for use in telescopes able to provide wide field imaging of 
the plasmasphere at 304 A. Utilizing a shutter to elliminate nonuniformities 
that occur during the interval when the mirror is being transported into and 
out of the deposition region will provide substantial further improvement. 
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Fig. 1 Photograph of the Multilayer Deposition System. 
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EXHIBIT A 


Effects of Doppler Shifts and Source Perspectives on Extreme Ultraviolet Images 
of Ion Populations Moving in the Inner Magnetosphere 

D. L. Murphy and Y. T. Chiu 

Lockheed Palo Alto Research Laboratory, Palo Alto, California 

Abstract. The full effects of solar EUV line resonant scat- 
tering by a population of ions in gyration, bounce and drift 
motion in the magnetosphere is formulated and demon- 
strated numerically with the 0 + 834 A line as test exam- 
ple. It is shown that anisotropies of motion induced by the 
geomagnetic field and source perspective effects caused by 
834 A earthshine profoundly affect the morphology and 
magnitude of the EUV line images, thus enhancing the 
effects of Doppler shifts on narrow EUV lines. The full 
effects, considered for the first time here, do not impact 
relatively wide resonant scattering lines such as the He + 

304 A line. 


Introduction 


The scientific concept and technology of imaging inner 
magnetospheric ion populations in solar extreme ultravio- 
let resonant scattering lines have been studied extensively 
in the last few years [e.g., Chakrabarti et al, 1982; Chiu 
et al, 1986, 1990; Meier , 1990; Swift et al. , 1989; see also 
the review by Williams et al, 1992 and references therein]. 
These independent estimates of signal intensity for some 
candidate resonant scattering lines were generally consis- 
tent with each other and yielded sufficient signal for mor- 
phological measurements, thus promoting the feasibility 
of the concept. While the above estimates are necessary 
to pioneer the concept, further definition of the charac- 
teristics is needed to examine how, and whether, the sci- 
entific objectives of imaging magnetospheric morphology 
and/or diagnostics of magnetospheric plasma processes can 
be achieved. This is the prime objective of this paper. 

This paper reports on results of imaging simulation 
studies aimed at complete inclusion of Doppler 
shifts with respect to illumination sources for resonant 
scattering by an ion population executing gyration, 
bounce, and drift motion in the geomagnetic field. As 
far as we can ascertain the effects of the geomagnetic field 
upon the EUV images of ion populations have not been 
treated in full, although Meier [1990] and Garndo et al 
[1991] have noted the importance of Doppler shifts on res- 
onant scatter rates. An important result of this study is 
that the anisotropic control of the geomagnetic field on ion 
population velocity distributions significantly impacts the 
morphology and diagnostics of some EUV images, espe- 
cially those with relatively narrow solar line widths. Fur- 
ther, sources of illumination strongly affect the perspective 
of some EUV images, because of the geomagnetic control 
' of ion motion. 
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Theory 

The theoretical treatment of resonant, scattering is de- 
scribed in detail by other authors [Mihalas, 1978; Meier, 
1991]. We briefly review the formalism to establish our 
notation. For magnetospheric ions we take the medium to 
be optically t hin. The radiance along a liue-of-sight ( LOS) 
r'(s) = r — sn is 


I{r, n i v)= f ds s(f\s), h,is) (L) 

d LOS 

where I is the specific intensity or radiance at observer 
position r for radiation of frequency v with direction of 
propagation n, e is the volume emissivity for such radi- 
ation, and 5 is the distance along the line-of-sight of the 
observer. There are three main ingredients necessary for 
evaluating the emissivity in resonant scattering: the source 
intensity, the amount of scattering for an individual ion, 
and the ion density. For a single resonant scattering line, 
the emissivity is 


c(r, h } is) = / — - dis'd 3 v J(r i ?, h, is, n , is') (2a) 

4x " J 4 7T 


with 


J(r, v , h, is , h\ is') = f(r, v)R(h y is , h\ is ' , F) J 0 (r, h\ is'). 

m 

In (2a) and (2b) /(r, F) is the phase space density of scat- 
terers, is 0 is the center frequency of the line, D \ 2 is the 
Einstein absorption coefficient, R is the Hummer redistri- 
bution function, normalized so that 


/ 


dQ, . dQ f _ , _ 
—dv—dis'R= 1 , 
4tT 4 7T 


( 3 ) 


and fo is the incident intensity. 

The formalism in this paper is applicable to any reso- 
nant scattering line, such as the helium 304 A line and the 
oxygen 834 A lines. Our results show that motion in the 
geomagnetic field does not affect the apparent morphology 
of 304 A images unless the helium ions are very hot, be- 
cause of the large solar line width; whereas it profoundly 
affects the oxygen 834 A images. Therefore, we formulate 
here the more complex oxygen problem in some detail. 

The solar 834 A complex comprises 9 spectral lines, a 
mixture of 0 + and 0 ++ emissions. We associate the val- 
ues of the index j from 1 to 3 with the three 0 + lines, and 
let values of j ranging from 4 to 9 correspond to the six 
0 ++ emissions. Meier [1990] has tabulated the statistical 
weights with which these lines appear in the solar 834 A 
blend. We use the notation isQj to denote the center fre- 
quency of the j th line, and Aisj to denote its width. We 
also denote by pj the probability that a solar photon in 
the 834 A line complex was emitted via the j th transition. 


15 


The source intensity is the sum over the three solar 0 + 
lines, the six solar 0 ++ lines, and the three earthshine 0 + 
lines. We first consider the solar emissions; we will discuss 
earthshine emissions below. For the solar lines we have 

I 0 = sf* hv F 0 fi{» ■ k - l) (4) 


In this expression, k is the direction of propagation of a 
solar photon and Fq = 6.7 x 10 8 photons/cm 2 /sec is the 
solar 834 A flux [Hintereggcr ct al , 1981]. We ignore the 
finite angular size of the sun. 

Following Meier [1991], Hummer’s redistribution func- 
tion for the k th magnetospheric 0 + line is given by 

H k = &[i>' - i/Qk(l + — — )] 6[v' - V + Vok -d- — — ] (5) 
c c 

if we assume an isotropic phase function. We will take 
the sum over the three 0 + 834 A lines below. The inci- 
dent radiation is described by frequency v ' and direction 
of propagation h f , the scattered radiation is described by 
the corresponding unprimed variables. 

The phase space density for magnetospheric 0 + is taken 
to be Maxwellian at an altitude of 1000 km, for all in- 
variant latitudes and magnetic local times. The model of 
Tsyganenko [1987] defines the magnetic field. We use con- 
servation of energy and first adiabatic invariant to map 
this distribution up magnetic field lines to other points in 
the magnetosphere; these two conservation laws will also 
provide the boundaries of integration in velocity space. For 
magnetospheric ion populations, the temperatures parallel 
and perpendicular to the magnetic field are generally not 
equal and the general formalism of Chiu and Schulz [1978] 
applies. If 7] | and T± are not equal, the resonant scatter- 
ing Doppler shift effects are more severely anisotropic with 
respect to the magnetic field direction. For the present 
work, we shall, for the sake of brevity of expressions, as- 
sume Tjj = 7T, which has the least anisotropy; we denote 
the common temperature by T. Thus, 


1 3/2 


/ = «0 


m 


2?r kT 


ex p[— 2]tr( v — v o)l 


( 6 ) 


where m is the mass of an 0 + ion, r 0 = |ro| = 1.1 57/Zg, 
and Vq is given by 


vl = 2GMe (---)- —Wi - *(*>)] (7) 

\r r 0 J m 

In (7) r 0 denotes the intercept of the magnetic field line 
with 1000 km altitude and <p is the magnetospheric electric 
potential. Although <j> is not utilized in the present work, 
which deals with a trapped 0+ ion population, we plan to 
use this capability in the future to incorporate the effects of 


parallel electric fields upon ion outflows. The boundaries 
of integration in velocity space are given by v 2 > r 2 and 
r 2 |fl(ro)| sin 2 or < [ v 2 — where a is the pitch 

angle, i.e., cos a = v B/(\v\\B\). 

With all ingredients at hand from (1), (5), and (b), the 
scattered intensity at the detector is given by substitu- 
tion of (2a-b) into (1). The intensity function /(r, n, v) 
is usually integrated over the instrument band pass. In 
this paper, we are not concerned with instrument effects; 
therefore, we integrate over the entire band assuming that 
there are no other background lines. 

The result for the total radiance, including the sums over 
spectral lines, is 


J 1 * 


1 2k 


/ 


where 


dsdPv J (1 + - — — ) exp [—A'] 
c 


v • k , 


A" = | ( ^0; ) /Ai/j 


(8a) 

( 86 ) 


Upon substitution for /, the integral over velocities be- 
comes, up to a constant of proportionality, 

J d 3 v (1 + l -^-) exp [—.4], (9) 

with 

A = Wr [v2 ~ v ° ] + x (10) - 

Using spherical coordinates i 1 , a, and / 3 in velocity space, 
the velocity space integration boundaries cited above can 
be expressed as step functions 

0[ir — vq] 0[(i> 2 — Vq)B — v 2 Bo sin 2 a] (11) 

with B — j#(r)| and Bo = |S(ro)|- The velocity space 
coordinate axes depend upon the direction of the magnetic 
field B(r). It is now clear that the integrated scattered 
intensity J I civ depends on direction with respect to the 
magnetic field upon velocity space integration. 

We include earthshine effects by assuming that above 
1000 km the illuminated hemisphere of the earth is a Lam- 
bertian radiator of 0 + 834 A emissions with a brightness of 
1.5 kilorayleigh. Note that we do not address the problem 
of the optically thick upper atmosphere, but only include 
its emissions as the earthshine source intensity illuminat- 
ing the magnetosphere above 1000 km altitude. The ra- 
diating atmospheric ions are assumed to have a thermal 
velocity distribution with a temperature of 1000 I\. This 
reduces the width from the solar width by the square root 
' of the ratio of the upper atmospheric temperature to the 
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solar atmospheric temperature. An additional complica- 
tion is the integration over the solid angle subtended by 
the illuminated earth, since the eartlishine source inten- 
sity is a non-trivial Function of the position along the line 
of sight and the photon direction of propagation. Thus, 
the evaluation of resonantly scattered eartlishine radiance 
involves a six dimensional integration-- the line of sight, 
the three dimensional velocity space, and the two angu- 
lar coordinates for the visible illuminated earth -whereas 
the resonantly scattered solar radiance is given by a four 
dimensional integration. 

Discussion and Conclusions 

We have evaluated these radiance integrals through 
Monte Carlo integration. In this integration method, the 
integrand is factored into two pieces. One piece is a prob- 
ability density function defined on the volume of integra- 
tion, e.g., the three-dimensional velocity space and the 
one-dimensional line-of-sight, for the case of solar illumina- 
tion. The other piece then becomes a function of random 
variables defined on the volume of integration, so that the 
original integral is now to be interpreted as the expectation 
value of this function. In Monte Carlo integration, a series 
of statistically independent points in the volume of integra- 
tion is chosen at random according to the aforementioned 
probability density function. The function whose expecta- 
tion value is to be determined is evaluated at each of these 
chosen points, and this ensemble of function values, when 
averaged, provides an estimate for the expectation value. 
By carrying out this process several times, an estimate for 
the uncertainty in the evaluation of the integral can also 
be obtained. For more information on Monte Carlo inte- 
gration, see Press et ai [1986] and references cited therein. 

We assumed a 10 eV oxygen ion reservoir at 1.157 i?£, 
with a spatial density of 100 0 + ions per cm 3 . The images 
shown are 31 x 31 pixels, smoothed to 60 x 80 pixels by 
linear interpolation. Statistical errors in the evaluation of 
the 0 + column density are about 1% or less for each pixel; 
statistical errors in the evaluated radiance range from 3% 
to 10%, with a typical value of 4%. We believe these er- 
rors to be uncorrelated from one pixel to another (prior 
to smoothing). With this integration code, we have simu- 
lated views in the magnetospheric 834 A band from 9 Re 
out for three spacecraft perspectives: tail ward, dawnward, 
and toward the direction of the negative ^-axis (in solar- 
ecliptic coordinates). Hereafter these will be referred to as 
‘dm ages' 1 for brevity. 

As has been pointed out by Meier [1990] and Gamdo 
et ai [1991], Doppler efFecls are important for 0+ 834 
A scattering in the magnetosphere. They appear in the 
resonant scattering formalism through the v • h and v ■ h l 
terms in the Hummer redistribution function. In order to 
highlight the effects of the Doppler shift, we compare the 
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full effects images with images produced with those terms 
set to zero. The image in Fig. la is a dawnward view, 
from 9 /?£, with all effects (Doppler shift, geomagnetic 
field, and earthshine) accounted for. Fig. lb shows the 
difference between an image including all effects except for 
Doppler effects, and the image in Fig. la. This difference is 
itself several times brighter than the full effects image, out 
to a few Re on the dayside. On each of these images, as 
well as on the images described below, we have marked on 
an outline of the earth the location of the magnetic north 
pole. 

We determine the importance of geomagnetic field ef- 
fects on magnetospheric images in EUV by evaluating the 
radiance integrals without the limitations set by geomag- 
neticaily trapped motion. The result of the angular inte- 
gration in velocity space is then scaled to give the same spa- 
tial density of 0 + ions, at each point along the line of sight, 
as would he obtained with the velocity space boundaries 
present. Thus, all pitch angles are available to the scat- 
tering 0+ ions, unlike the geomagnetically trapped popu- 
lation. To our knowledge, the evaluation of the influence 
of geomagnetic effects on EUV images has not heretofore 
been carried out. The basic cause of the geomagnetic ef- 
fects is the anisotropy induced by the geomagnetic field 
on the parallel and perpendicular velocity distributions of 
a trapped population executing gyration and bounce mo- 
tion. With all pitch angles available, the expected radi- 
ance is enhanced in the auroral regions and reduced in the 
equatorial regions. Unless the importance of geomagnetic 
field effects on 0 + ion pitch angle distributions is recog- 
nized, misinterpretation of such magnetospheric EUV im- 
ages will result. Figs. 2a and 2b display the differences in 
morphology due to the inclusion/omission of pitch angle 
effects. Fig. 2a is an image of the difference: [ail effects] 
- [all effects but for pitch angle effects]; Fig. 2b displays 
the difference: [all effects but for pitch angle effects] - [all 
effects]. In both images, values of the radiance difference 
less than 10 -3 Rayleigh are set to white. Comparing these 
difference images with the full image Fig. la, we see that 
(1) the magnitude of the difference is a substantial part of 
the full image at the important altitudes > li?£, and (2) 
the difference images in Fig. 2b have the morphology of 
magnetic field lobes, showing the important combined in- 
fluences of Doppler shifts and geomagnetic field on imaging 
the morphology of the 0 + population. 

Earthshine is an important source of 834 A illumination 
out to a couple Re on the dayside, as shown by Fig. 3b, 
which is to be compared to Fig. 3a. Each of these im- 
ages is a view looking toward the negative solar-ecliptic z 
axis from a distance of 9 Re- Fig. 3a includes all effects, 
whereas Fig. 3b omits the effects of earthshine and includes 
only Doppler and geomagnetic effects. The relative impor- 
tance of solar and terrestrial sources of EUV depends upon 
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solar and magnetospheric activity and atmospheric condi- 
lions. Comparison of Figs. 3a and 3b shows that, whereas 
the morphologies of the two are similar, the Doppler and 
geomagnetic effects combined with earthshine effects pro- 
foundly influence the expected image intensity, thus im- 
pacting the proper interpretation in terms of plasma diag- 
nostics. 
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Figure Captions 


Fig. I: (a) Simulated magnetospheric 854 A image from 
9 Rt ; above the dusk meridian, looking dawnward, includ- 
ing effects of Doppler shifts, velocity space anisotropies due 
to the geomagnetic field, and earthshine; (b) same perspec- 
tive as in (a), but now the image shown is the difference 
between an image omitting Doppler effects and the image 
in (a). 


Fig. 2: Views displaying the difference between a simulated 
image including all effects, and one omitting geomagnetic 
pitch angle effects, for the perspective of Fig. 1. (a) [all 
effects included] - [geomagnetic effects omitted]; (b) [geo- 
magnetic effects omitted] - [all effects included]. 


Fig. 3: Simulated magnetospheric images looking down 
from the positive solar-ecliptic z-axis to show the effects 
of earthshine. (a) is an image with all effects included; in 
(b), the effects of earthshine have been omitted. 
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EXHIBIT B 


Magnetospheric Imaging of High Latitude Ion Outflows 

D.E. Garrido', R.M. Robinson, Y.T. Chiu, H.L. Collin 

Lockheed Palo Alto Research Laboratory, 3231 Hanover St., Palo Alto, CA., 94304, U S A. 

R.W. Smith and D.W. Swift 

Geophysical Institute and Department of Physics, University of Ala ska- Fairbanks, Fairbanks, AK . 99775-0800, U.S.A 

Simulation of magnetospheric emissions from resonant scattering of sunlight by upward flow- 
ing /fe + (30.4 nm) and 0 + (83.4 nm) ion» at high latitudes have been constructed for different 
satellite viewing locations and distances from the earth. The ion outflow data obtained by the 
Lockheed Ion Mass Spectrometer (EICS) on board the DEM spacecraft has been used to take 
Into account the non-uniform distribution of ion flux intensities. For a polar orbiting satellite we 
obtain emission rates from outflowing ions of 5.2xt0~ a and 5.flxI0~* Rayleighs for quiet and 
disturbed times for He + and 8.5x10” 3 and 7.5xlO~ J Rayleighs for quiet and disturbed times 
for 0 + , respectively. These low light levels are due to use of ion outflow data only. Emission 
levels from all iom can be quite a bit higher and awaits further data reduction from the Lock- 
heed instrument. We find that viewing at a distance of 9 Hjj with an instrument field of view 
of at least 30 degrees that includes an occultation disk in the image is sufficient to exclude the 
bright emissions coming from plasmaspheric He + ions. Because the 0 + intensities increase dra. 
matically with magnetic activity, they represent a good signature of subatorm-related processes. 

In our modeling we assume that the ions move along dipole magnetic field lines with a density 
profile that is modeled by a simple magnetic field scaling. The Doppler shift of these moving ions 
modify the intensity of the resonantly scattered emissions. For //e + the Doppler effects tend to 
decrease the intensities of the magnetospheric images while for 0 + at sufficiently high velocities 
the intensities may increase because enhanced scattering occurs. 

1. INTRODUCTION 

Interest in global imaging of the magnetosphere in the ul- 
traviolet has been initiated by Swift et al. (1989), Chiu et 
ah (1990), Meier (1991) and Garrido et al. (1991). Swift 
et al. (1989) constructed images of the plasmasheet and 
ring currents from resonantly scattered light from 0 + ions. 

Chiu et al. (1990) simulated images of outflowing G + (83.4 
nm) ions and neutral O (130.4 nm) atoms by considering 
a geostationary satellite viewing at a distance of 6 R& in 
the dusk meridian. They assumed particle motions along 
a static dipole field with tilt and a density profile that is 
modeled by a simple magnetic field scaling. A uniform flux 
distribution over the polar region was assumed with a value 
derived by averaging outflowing ion fluxes measured by the 
EICS instrument on DE-1. The resulting column emission 
rates ranged between .01-1 Rayleighs. Effects of gyromo- 
tion of the ions and Doppler shifting due to the ions bulk 
velocity along the sun-earth direction was not considered. 

The calculations of Chiu et al. (1990) also did not include 
the bright resonant emissions originating from the plasma- 
sphere which can obscure those emissions coming from the 
high latitude ion fluxes. 

Meier (1991) provided a detailed assessment on the pos- 
sibilities of remote sensing in the ultraviolet over a broad 
range of wavelengths that includes the 30.4 and 83.4 nm 
resonant lines of He + and O*, respectively. Using plasmas- 
pheric models based on the observations of Meier and Weller 
(1972) and Weller and Meier (1974), images of the plasma- 
sphere were constructed at a viewing distance of 10 Re in 
the dawn to dusk meridian plane and at various satellite 
locations defined by geographic latitudes 0, 30, 60 and 90 
degrees. The extent of the plasmapause boundary is visible 
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A. PROGRESS 


This report summarizes work performed to date on NASW-4529: Imaging the 
Magnetosphere in the Extreme Ultraviolet. In the ten month's elapse time, contract activity 
was focused on three tasks scheduled to proceed more or less in parallel. These are: 

(1) Establishment of reliable estimates of expected solar EUV resonance scattering fluxes 
at He + 304A and 0+834A lines based on He + and 0 + fluxes measured by Lockheed s ion 
mass spectrometer (EICS) on board the NASA Dynamics Explorer satellite. 

Mr. Harry Collin, ion mass spectrometer data scientist in our laboratory, participated in this 
data reduction, effort as did Mr. Dante Gamdo, graduate student from the University of 
Alaska Fairbanks working on magnetosphere imaging as thesis topic. Under agreement by 
the PI and Dr Akasofu of UAF, Mr. Gamdo was sponsored by Lockheed via UAF s 
NASA Space Grant program for which Lockheed agreed to provide salary support to Mr. 
Garrido to participate in the data analysis and the computer image simulation efforts. 
Because this data analysis effort requires more resources than is allocated in NASW-4529, 
part of the data analysis effort is funded by Lockheed's independent research program. 

This collaborative effort proves to be fruitful as Mr. Garrido's work at Lockheed will 
constitute the major part of his Ph.D. thesis and it leads to a publication. 

Two publications resulted from our effort to simulate magnetosphere EUV imaging based 
on in situ ion mass spectrometer data. The concept is to map all relevant in situ ion mass 
spectrometer data along magnetic field lines to an altitude of 1000 km in order to build up a 
statistical global distribution of ions on the 1000 km surface. Then the statistical fluxes 
, . distribution is treated as an instantaneous ion flux distribution, which we map up the field 

line as a surrogate magneto spheric ion population. Computer simulation of a 
magnetosphere image then proceeds by allowing the ions to resonantly scatter solar EUV 
lines into an imager of given optical characteristics. The results are summarized m the 
abstracts of the two publications whose front pages are attached to this report as Exhibits 1 
and 2. Exhibit 1 is published in the Geophysical Research Letters and Exhibit 2 is to be 
published in Planetary and Space Science, in deference to our UAF colleagues' preference. 

(2) Multilayer mirror design for a wide field-of-view 304A imaging telescope. 

For space-based imaging the magnetosphere in the next decade, our examinations, with 
computer simulations such as those in Exhibits 1 and 2, indicate that the only viable 
technique is that of multilayer mirror normal incidence imaging in the He + 304A line. No 
834A multilayer mirror has been manufactured or flown to date. The first and only 304A 
mirror has been flown by a Lockheed-Stanford collaboration for astrophysical objects 
imaging, although such imaging requires relatively small field-of-view. For wide field-of- 
view imagery required for magnetosphere imaging, multilayer mirror design has to be 
considered afresh. Dr. Dick Catura of our laboratory is the leader of this multilayer mirror 
design and manufacture effort for the astrophysical collaboration. 

Dr. Catura was asked to consider design of an imaging mirror at 304A with large field of 
view. Exhibit 3, the front page of a paper by Dr. Catura, illustrates the capability and 
facility of our laboratory on multilayer mirror design and manufacture. 

A molybdenum-silicon multilayer mirror suitable for 304A magnetosphere imaging has 
been designed with good angular response of ± 10 degrees. Exhibit 4 shows the reflectance 
^ ^ curve of the designed multilayer with a sharp peak of about 15% at 300A . Filters are 

usually required to cut out unwanted reflectance at both long and short ends of the 
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spectrum. Such filters can be designed as multilayers also. Exhibits 5 and 6 show the 
angular response of the multilayer for two options of slightly different multilayer 
parameters. The important point for magnetospheric imaging is that a wide field of view 
mirror is possible. 

(3) Extraction of physical parameters from images. 

For magnetospheric imaging to be a science project of sufficient value, it is necessary to go 
beyond the mere making of a magnetospheric image. Consequently, we devote a small 
portion of our resources this year to investigate the theoretical possibility of extracting 
information from the images. Among these, we have chosen to investigate the possibility 
of extracting ion temperature and/or ion outflow information in the inner magnetosphere. 
Further, to make reliable images, the Doppler shift of moving magnetospheric ions relative 
to the resonant scattering line frequencies must be investigated. This important point has 
recently been pointed out by our collaborator, Dr. Bob Meier of NRL, and a linear motion 
treatment was given by Dr. Meier for the case of the 0 + 834A line (Exhibit 7) . However, 
for a hot ion population gyrating and outflowing along a magnetic flux tube, the Doppler 
shift is not so simply accounted for. As a result, we initiated a proper calculation of 304A 
and 834A resonance scattering by He + and 0 + ions with proper magnetospheric 
characteristics. The calculation is quite complex, as can be seen from Exhibit 8, which 
summarizes its basic ingredients. 

B. FUTURE WORK 

For continuation in 1992 of the current project, we intend to proceed with the second phase 
of all three tasks whose progress in 1991 are summarized above. 

(1) Imaging simulation with Lockheed data. 

For this task, the 1992 emphasis will be to push the consideration of our publication 
(Exhibit 2) towards lower latitudes. This will be consonant with objectives of NASA’s 
Inner Magnetosphere Imager mission. The state of data processing for Lockheed s EICS 
instrument for lower latitudes is not as advanced as for higher latitudes. However, we 
estimate that sufficient resources will be available to complete this task in CY 1992-1993. 

(2) Multilayer mirror design for 304A imaging. 

The optical properties of our mirror designed will further be defined with computer 
simulations. A small mirror will likely be available from our parallel Lockheed independent 
research project which is not confined to space-based imaging. We intend to perform 
vacuum exposure and long-term outgassing experiments on this mirror in order to 
determine its performance in the space environment Finally, if resources are available we 
intend to revisit the critical issue of minor uniformity because it is critical to the upgrade 
from small field-of-view astrophysical imaging applications to wide field-of-view 
magnetospheric imaging applications. 

(3) Image Exploitation 

We shall complete the formulation of the Dopper shift problem outlined in Exhibit 8 of the 
PROGRESS section. A computer simulation of the problem outlined in Exhibit 8 will be 
completed in CY 1992. 
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GEOPHYSICAL RESEARCH LETTERS, VOL. 17, 


NO. 3, PAGES 267-270, MARCH 1990 


MAGNETOSPHERIC AND EXOSPHERIC IMAGING IN THE EXTREME ULTRAVIOLET 

Y. T. Chiu, R. M. Robinson, H. L. Collin 
Lockheed Palo Alto Research Laboratory, Palo Alto, California 
S. Ghakrabarti, G. R. Gladstone 
Space Sciences Laboratory, University of California, Berkeley 


Abstract. Extreme ultraviolet line emissions by exo- 
teric oxvgen ions and neutrals of ionospheric origin can 
3 e used to form images of the dynamic magnetospheric re- 
dons with sufficient speed to provide a new global means 
if observing the dynamics of magnetosphere-ionosphere 
•oupling processes. Computer simulations of such ex- 
jspheric images based on the latest measurements of 
on outflow, on known or estimated solar fluxes, and on 
mown emission and excitation rates are shown. An evab 
nation of the current and foreseeable EUV imaging tech- 
nology relevant to magnetospheric imaging is also given. 

Introduction 

The idea of using photometric or imaging techniques 
on ionic emission lines to obtain images or density dis- 
tributions of magnetospheric and exospheric regions of 
the Earth began with observations of far ultraviolet emis- 
sions from neutral hydrogen and ionized helium reported 
bv Young et al. 11971] and Johnson et al. [19 <11. 
Chakrabarti et al. |1982] used solar resonance scattered 
emissions from Hell at 304Ato map the helium distribu- 
tion in the plasmasphere. The full concept of EU 
ospheric imaging with a set of EUV lines was discussed 
together with a search for possible verifications from some 
existing photometric data by Chiu et al. [1986 ; instru- 
ment concepts were also discussed [Chakrabarti et al., 
19861 These studies indicated that exosphenc imaging 
was technologically feasible in the inner magnetosphere 
f<10 Rf) with what was known about emission sources 
and instrument capabilities at the time Recently, there 
has been a resurgent interest in EUV magnetosphenc 
imaging at greater distances using observations from a 

lunar base [Swift et al., 1989]. , 

The purpose of this paper is to re-examine the exo- 
spheric imaging concept in the light of ^ hree 

new developments: (1) new information about the outflow 
of ionospheric ions, (2) recent calculations of the escape 
fluxes of neutral 0° and H° due to charge exchange, and 
(3) the development of layered synthetic microstructure 
(LSM or multilayer) sensors which will result in drastic 
reductions in the size, increases in the speed, and enhance- 
ments in imaging quality of EUV instruments. 

M agnetospheri c and Exospheric Imaging 

Magnetospheric imaging requires the observation of 
emissions of ionized constituents because the distnbution 
of these particles is constrained by the magnetic held. 
During quiet times H + is the dominant ion in the magne- 

Copyright 1990 by the American Geophysical Union. 

Paper number 90GL00323. 

0094-8276/90/90GL-00323$03 . 00 


tosphere, but H + emissions cannot be detected optically. 
He^ is also a candidate for use in magnetosphenc imaging 
but there is not enough He + in the outer magnetosphere 
to image the regions of interest. 0 of ionospheric origin 
can have sufficient concentration in the magnetosphere for 
use in imaging. Also, escaping 0 + charge exchanges and 
produces 0° which can be imaged. Although the 1304.4 
triplet is optically thick in the lower thermosphere, at al- 
titudes above ~1000 km it is optically thin. Since the 
ground state sub-levels of the exosphenc 0 are not ex- 
pected to be in thermal equilibrium, only the 1302A com- 
ponent has been used in calculating the g-factor for this 

feature. , „ . , r 

The EUV line emissions that we shall consider for imag- 
ing are produced bv resonance scattering of solar light by 
the relevant constituents. In the following, si mu a lions 
and observational techniques for EUV images of 0 nul 
O 0 will be discussed as an example of the key elements ot 
magnetospheric imaging. Lack of space precludes discus- 
sion of other emission lines. . . . r 

The brightness, 4 t r/, of an optically- thin emission line 

excited by solar scattering is given by 


4tt7 


= <7 J N(t) 


dl 


( 1 ) 


where J N dl is the column density of atoms/ions along 
the line of sight l and the g-factor is estimated by 


g = TT 1 ^ 2 <T 0 A Xd * F\ 

~ <T 0 (AA$/AAq) 7rF 


( 2 ) 


Here <r 0 is the resonance-scattering cross section at line 
center A 0 , AA D is the Doppler width of the scattering 
atoms/ions defined by AA D = AqU^/c, , is the so- 
lar flux per unit wavelength at line center and rrl is the 
total solar flux of the line. For emissions such as the 
834A line, for which the solar line profile is unknown 
an estimate was made using a value of the Earth and 
solar Doppler width ratio AAq/AAq ~ 0.05. Table 
lists the total solar fluxes and estimated g- factors for line 
emissions. Cross sections used in Table 1 were cMculated 
from oscillator strengths found in Wiese et al. [1966]; the 
solar fluxes were estimated from Heroux and Hintereg- 
ecr 119781, and line shapes were estimated from bkellon 
f]()78|. Recent solar models [D. Sliemansky, private com- 
munications, 1989] suggest that perhaps most of the solar 
emission at 834A results from a coincident 0 ++ multiplet 
rather than from the 0+ 834A triplet. The estimated g- 
factor for the 834A line reflects this new information If 
this suggestion proves not to be true, then the 834 A g- 
factor in Table 1 may have to be increased by perhaps a 

Simple models of 0+ and 0° number densities have 
been used, based on observations and/or theory, with the 
contribution due to ‘hot’ outflowing ions and atoms ex- 
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Title: Magnetospheric Imaging of High Latitude Ion Outflows 

D.E. Garrido, R.M. Robinson, Y.T. Chiu, H.L. Collin 
Lockheed Palo Alto Research Labortories, Palo Alto, CA 

R.W. Smith and D.W. Swift 
University of Alaska, Fairbanks, Alaska 
Abstract: 

Simulation of magnetospheric emissions from resonant scattering of sunlight by upward flowing 
He + (30.4 nm) and 0 + (83.4 nm) ions at high latitudes have been constructed for different satellite 
viewing locations and distances from the earth. The ion outflow data obtained by the Lockheed 
Ion Mass Spectrometer (EICS) on board the DE-1 spacecraft has been used to take into account the 
non-uniform distribution of ion flux intensities. For a polar orbiting satellite we obtain emission 
rates from outflowing ions of 5.2 x 10" 2 R and 5.6 x 10' 2 R for quiet and disturbed times for He + 
and 6.5 x 10‘ 2 R and 7.5 x 10* 2 R for quiet and disturbed times for 0+ respectively. These low 
light levels are due to use of ion outflow data only. Emission levels from all ions can be quite a bit 
higher and awaits further data reduction from the Lockheed instrument. We find that an instrument 
field of view of at least 30° is required at a distance of 9 Re for images of the outflowing ions that 
exclude the bright emissions coming from plasmaspheric He + ions. 

Because the 0 + intensities increase dramatically with magnetic activity, they represent a good 
signature of substorm-related processes. In our calculations we assume that the ions move along 
dipole magnetic field lines with a density profile that is modeled by a simple magnetic field scaling. 
The doppler shift of these moving ions modify the intensity of the resonantly scattered emissions. 
For He + the Doppler effects tend to decrease the intensities of the magnetospheric images while for 
0+ at sufficiently high velocities the intensities may increase because enhanced scattering occurs. 
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Fabrication and Analysis of Mo-Si Multilayers 
L. Shing, R.C. Catura 
0/91-30, Bldg 256 

Lockheed Palo Alto Research Laboratory 
3251 Hanover St. 

Palo Alto, CA 94304 


ABSTRACT 

Molybdenum- Silicon (Mo-Si) multilayers have been fabricated in our laboratory for reflecting various 
soft x-ray wavelengths. Several normal incidence mirrors have been produced for a solar sounding rocket 
flight. Using properly shaped masks on 3 inch diameter DC magnetrons, we can produce multilayers of 
a similar size having a uniformity of better than ±3% over their entire surface area. Peak reflect^ce of 
the multilayers are as high as 42% at 171 A. Efforts are underway to deposit several layers of Afferent 
period on top of the multilayers in order to attenuate the reflection of 304A, a very bright line of He II 
in the solar Spectrum. In the interim, we have made a Mo-Si mirror that reflects at 304A. In addition 
to obtaining information on the solar corona at this wavelength, this will greatly reduce the effects of the 
He II line when its image is appropriately subtracted from others. These results and a discussion of he 
current status of our multilayer research will be presented. 

2 INTRODUCTION 

We have been involved in fabricating multilayers for several years and our most recent experience is in 
depositing Molybdenum-Silicon multilayers on figured optical surfaces. First, the performance of various 
multilayer designs are studied by computer simulation to determine the optimal design ^ a ^elected 
wavelength, bandpass or angle of incidence. Step heights of each material are then deposited onto standard 
Si wafers to calibrate the deposition rates of the magnetrons. Proto-type multilayers are deposited onto 
Si wafers to confirm integrity of the design, and adequacy of the fabrication setup. The reflectance o e 
proto-types are then measured at various wavelengths, angles of incidence, and position on the test wafer. 
After optimizing the physical parameters, the actual flight optic is coated, and if practical, its reflectance 
U also measured. This paper describes our most recent experience in fabncatmg Mo-Si multilayer optics 
and presents data on their measured performance. 

3. COMPUTER SIMULATION 

We are able to calculate the reflectivities of various multilayer designs using a computer code developed 
by Dr David Windt and his compilation of optical constants for a large number of material. The code 
is able to compute the response of both periodic and nonperiodic multilayers. Figures 1 and 2 show the 
calculated reflectivities of multilayers designed for use in solar observations during an upcoming sounding 
rocket flight 1 . Parameters of the multilayer designs are indicated on the plots. The multilayer passband 
at 17lA is of interest to observe structure in the solax corona in an emission line of ionized Fe (Fe IX) tha 
is formed at temperatures near 10 a K. That, at 195 A w.U be used to observe structure in a hue of Fc XII 

formed at coronal temperatures near 2 x 10 K. , . a • 

The computer simulations of the response of these multilayer designs show that their reflectivities 

at 304A are about 5% of the reflectivities at the wavelengths of interest. This attenuation, however, is 
insufficient to reduce the contribution of an intense solar emission of He 11 at 3°4A, which is about two 
orders of magnitude brighter than the two emission lines of interest. A possibility for further attenuating 
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